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ABSTRACT 

We test the gamma-ray burst correlation between temporal variability and peak luminosity of 
the 7-ray profile on a homogeneous sample of 36 Swift/BAT GRBs with firm redshift determi- 
nation. This is the first time that this correlation can be tested on a homogeneous data sample. 
The correlation is confirmed, as long as the 6 GRBs with low luminosity (< 5 x 10 50 erg s _1 
in the rest-frame 100-1000 keV energy band) are ignored. We confirm that the considerable 
scatter of the correlation already known is not due to the combination of data from different 
instruments with different energy bands, but it is intrinsic to the correlation itself. Thanks 
to the unprecedented sensitivity of Swift/BAT, the variability/peak luminosity correlation is 
tested on low-luminosity GRBs. Our results show that these GRBs are definite outliers. 

Key words: gamma-rays: bursts - methods: data analysis 



1 INTRODUCTION 

A number of correlations between intrinsic properties of Gamma- 
Ray Bursts (GRBs) has been discovered since it has become pos- 
sible to measure their distances. In particular, correlations between 
properties of the 7-ray prompt emission as well as of the afterglow 
at different wavelengths have provided an increasing number of 
clues to identify the mechanisms and, ultimately, the nature of the 
GRB progenitors. In addition, some of these correlations have been 
tentatively used as luminosity estimators, with several implications 
on their possible usage to constrain the cosmology of the Univers e 
dGhirlanda et alj2004l ; iLiang & Zhangll2005l ; iFirmani et af]|2005h . 

The increasing number of GRBs with spectroscopic redshift 
allows to test and better calibrate them. Rece ntly, a crucial contri - 
bution has been supplied by the Swift satellite (Geh rels et alj2004t) . 
whose average rate of 100 GRBs per year since launch (November 
2004) made it possible to measure the distances of almost 1/3 of its 
sample, thus duplicating the overall number of GRBs with known 
redshift since 1997. 

The sample of GRBs detected with the Swift Burst Alert Tele- 
scope (BAT; Barthelmy et al. 2005) is particularly suitable to test 
the correlations between intrinsic properties, with the unprece- 
dented benefit of a homogeneous data set, apart from those requir- 
ing the peak energy measurement, made difficult by the limited en- 
ergy band (15-350 keV). 
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Hereafter we focus on a long-standing correlation between 
the variability and peak luminosi t y of the 7-ray promp t emission 
dFenimore & Ramirez-Ruizj 2000; Rei chart et alfe oOl). In partic- 
ular, Reichart et al. (2001; hereafter R01) provided a definition of 
variability (hereafter denoted as Vr) that turned out to correlate 
with the isotropic-equivalent rest-frame 100-1000 keV peak lumi- 
nosity (hereafter L) for a sample of 1 1 GRBs with known redshift 
availa ble at the time, usin g data from the CGROfBKTSE experi- 
ment jPaciesas. et alJl999T) . 

R01 modelled the variability/peak luminosity correlation 
(hereafter V/L correlation) with a power law (L oc V^™) with 
m — affected by extrinsic or sample scatter, described 

by criogVp, = 0.18. Recently, Guidorzi et al. (2005; hereafter 
GFM05) and Guidorzi (2005; hereafter G05) tested the V/L cor- 
relation on an extended sample of 32 GRBs with known redshift 
(GFM05) and on 551 BATSE GRBs, respectively. For the latters, a 
pseu do-redshift was derived assuming th e lag-luminosity correla- 
tion jNorris et al.ll2000l : lBand et all2004h . 

Both works confirmed the correlation, but with a lower slope 
than that derived by R01: m = 1.3±%* (GFM05) and m = 
0.85 ± 0.02 (G05). However, in either case it was pointed out 
that the scatter around these pow er laws made the description of 
a simple power law unsatisfactory. iReichart & NvsewandeJ J2005h 
applied the same method as R01 to the very results obtained by 
GFM05, obtaining m = 3.4± ;^ and ai og v B = 0.20 ± 0.04, per- 
fectly in agreement with the original values of R01. They ascribed 
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the disagreement to the fact that GFM05 did not deal with the sam- 
ple variance properly. 

More recently, Guidor zi et alj d2006h applied the D'Agostini 
(2005) method, accounting for the sample variance, to the data 
sets of both GFM 05 and G05. They obtained shall ower slopes than 
those by R01 and lReichart & Nvsewandeil |2005) and larger scat- 
ters: in particular, for the sample of 32 GRBs with firm redshift 
drawn from GFM05 they obtained m = 1.7±0.4, <nogv R ~ 0.34, 
while for the sample of 551 GRBs with pseudo-redshifts of G05 it 
resulted m = 0.88±£", cri og Vr ~ 0.74. 

For more details on the debate concerning the methods 
to be used, we refer the r eade r to the original p apers by 
iReichart & Nvsewandeil d2005h and lGuidorzi et alj d2006l) . 

Li & Paczyhski (2006; hereafter LP06) have recently provided 
a slightly modified definition of variability, hereafter denoted as 
VLp, which they found to correlate more tightly with L than Vr, 
without any extrinsic scatter in addition to the uncertainties affect- 
ing the single values of the single GRBs. Vlp differs from Vr 
mainly in the choice of the smoothing filter determining the ref- 
erence light curve with respect to which the variance is evaluated. 
LP06 chose the Savitzy-Golay filter instead of a simple boxcar used 
by R01. As a result, Vlp selects only the high frequencies, whereas 
only in the Vr calculation the lower frequency variance can give a 
contribution. 

The variability of the 7-ray prompt emission light curves is 
supposed to be produced above the photospheric radius of the fire- 
ball, above which radiation becomes optically thin. The interpre- 
tations proposed of the V/L correlation mainly invoke the pres- 
ence of a jet, whose angle 0, i.e. either the opening angle or the 
viewing angle (e.g., see Ioka & Nakamura 2001) for some jet pat- 
terns, is strongly connected with the observed peak luminosity L 
as well as with the Lorentz factor V of the expanding shell(s). The 
result would be a strong dependence of both L(9) and T{6) on 8. 
For instance. iKobavashi et al.l J2002h reproduced the observed cor- 
relation through numerical simulations, assuming F cx 9~ q and 
a log-uniform distribution in the time delay between next shells, 
from 1 ms to 1 s. A value of q = 2 seems to account well for 
the results bv lGuidorzietal.ld2006t) as well as the anti-correlation 
betwe en break time and peak luminosity dSalmonson & Galamal 
|2002|). Similar results have been found by iMeszaros e t al. ( 2002) 
and Ramirez-Ruiz & Llovd-Ronnind d2002l) under slightly differ- 
ent assumptions. 

The new piece of information from this analysis is given by 
the presence of low-luminosity high-variability GRBs. 

In this paper, we test the V/L correlation on a homogeneous 
sample of 41 GRBs detected with Swift/B AT using fully homoge- 
neous data. We considered two different definitions of variability: 
that by R01 and that by LP06. In Section [2] we describe the data 
sample and the selections we made. Sections [3] and [4] report how 
peak luminosity and variability have been calculated. Results are 
reported in Sec.[5]and discussed in Sec. [6] 



2 THE GRB SAMPLE 

The sample includes 5 1 long (T9 > 2s) GRBs with spectroscopic 
redshift detected by Swift/B AT dGehrels et al]|2004|) between the 
launch (2004, November 20) and October 2006. Out of this sam- 
ple we selected only those bursts who se 7-ray profile is ent irely 
covered by BAT during the burst mode dBarthelmv et alj2005h . No 
further selection was made on the sample, in order to avoid any ar- 
bitrary bias in the results. This requirement resulted in the rejection 



of 10 GRBs. In fact, in these cases the observation of Swift/BAT 
switched from burst mode to the survey mode before of the end of 
the prompt emission. The light curve results with a truncated pro- 
file. This is the case of GRB 050318, whose light curve stops about 
32 s after the trigger, as well as of GRB 050820A, GRB 050904 and 
GRB 060218. For GRB 060124 only the precursor was recorded in 
event mode, while the main event was observed in survey mode. 
For GRB 060906 the light curve is incomplete at the beginning, 
because the trigger probably missed the true onset of the burst. 
No burst mode event file is available for GRB 060505, as BAT 
observed it only in survey mode. We chose not to make use of 
the background subtracted light curves acquired during the survey 
mode to keep the sample as homogeneous as possible. GRB 050408 
was detected by XRT and UVOT, but not by BAT, although the 
light curve of its prompt emission is available from other instru- 
ments {HETE-2/FREGATE; Atteia et al. 2003). Nevertheless, we 
did not consider it in this work because we focused on BAT data 
for the reasons reported above. In the case of GRB 050802 and 
GRB 051227A the problem is in the red shift determination. Fo r 
the former only a tentative redshift exists dCummings et alj|2005l) . 
which is at odds with t he interpretation of the SwiftfUVCTT results 
dMcGowan et alj|2005h . For GRB 0512 27A there is a reds hift de- 
termination of the putative host galaxy dFolev et al]|2005j) . but it 
is still unclear if this is the real host galaxy. 

After this selection the sample has shrunk to 41 long GRBs, 
entirely covered by BAT and processed through the same proce- 
dure. Therefore, this work investigates the V/L relation based on a 
completely homogeneous sample. 

The BAT event files were retrieved from the Swift public 
archive and analysed through the standard BAT analysis soft- 
ware distributed within FTOOLS v6. 1. For each GRB we extracted 
mask-tagged light curves for a number of different binning times 
in the total nominal energy band (15-350 keV) 0, through the 
tool batmaskwtevt adopting the ground-refined coordinates pro- 
vided by the BAT team for each burst. These curves are there- 
fore already background subtracted according to the coded mask 
technique (Barthelmy et al. 2005 and references therein). For each 
burst the BAT detector quality map was obtained by processing the 
next earlier enable/disable map of the detectors, telling which de- 
tectors were disabled in flight because too noisy. We also applied 
the energy calibration to the event file making use of the closest-in- 
time gain/offset file through the tool bateconvert, as suggested by 
the BAT tearrQ. Finally these light curves are expressed as count 
rates with uncertainties: the rates are background- subtracted counts 
per second per fully illuminated detector for an equivalent on-axis 
source, as the default corrections are applied: ndets, pcode, maskwt, 
flatfield. 

We also studied the behaviour of the background fluctuations 
in burstless regions of the light curves and we found that the mask- 
tagged rates, r\, fluctuate compatibly with a white noise with sigma 
a Ti (ri and a Ti are the rate and its uncertainty of the i-th bin, 
respectively; see Appendix). We concluded that an upper limit of 
~ 2-4% (4-6%) at 90% (99%) confidence level can be derived on 
the presence of a possible extra variance (of instrumental origin, for 
instance) in addition to that due to the Poisson counting statistics, 



1 http://swift.gsfc.nasa.gov/docs/swift/archive/ 

2 The effective band is 15-150 keV, because photons with energy above 
150 keV become transparent to the coded mask and are treated as back- 
ground by the mask-weighting technique (e.g., Sakamoto et al. 2006). 

3 http://swift.gsfc.nasa.gov/docs/swift/analysis/threads. 
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implicitly assumed during the light curve extraction with the tool 
batbinevt. 

We found that it is not correct to perform the same analy- 
sis on BAT light curves with raw counts, i.e. not masked. In fact, 
we found that the GRB profile itself can be dramatically contam- 
inated by other sources and by background variations, with time, 
due to the slewing of the spacecraft during the prompt emission, for 
most GRBs. Furthermore, we found that BAT light curves with raw 
counts are severely affected by extra variance, which is comparable 
with the Poisson variance due to the counting statistics, in agree- 
ment with previous results (LP06). Therefore, we conclude that the 
BAT light curves of most GRBs with raw counts, not masked, are 
not suitable for temporal variability studies. 



3 PEAK LUMINOSITY 

For each GRB we extracted the mask-tagged light curve with a bin- 
ning time of 50 ms in the 15-350 keV energy band. We determined 
the 1-s time interval with the highest total counts and assumed this 
as the time interval corresponding to the 1-s peak count rate. 

We extracted the mask weighted spectrum in this time interval 
using the tool batbinevt. We applied all the corrections required: 
we updated it through batupdatephakw and generated the detector 
response matrix using batdrmgen. Then we used batphasyserr in 
order to account for the BAT systematics as a function of energy. 
Finally we grouped the energy channels of the spectrum by impos- 
ing a 5-cr (or 3-<7 when the S/N was too low) threshold on each 
grouped channel. We fitted the resulting photon spectrum, $(-E) 
(ph cm _2 s _1 keV _1 ), with a power law with pegged normalisation 
(PEGPWRLW model under XSPEC v. 12), except for GRB 050525A 
and GRB 060927 where a cutoff power law was used, in the rest- 
frame energy band 100 — 1000 keV. The choice of the energy band 
is connected with the original definition by R01 (see their eq. 9) 
also used by GFM05 (their eq. 7). 

Therefore the GRBs rest-frame 100 - 1000-keV isotropic- 
equivalent peak luminosities were computed using: 



(.1000/(l+z) 

4ttD 2 l (z) / E<$>{E)dE 

Jl00/(l+z) 



(1) 



where Dl(z) is the luminosity distance at redshift z, E is energy 
expressed in keV. Finally we derived the uncertainty on the peak 
luminosity by propagating that of the measured flux. 

Concerning the six BAT GRBs shared with the sample of 
GFM05, we compared the two sets of peak luminosities: these 
GRBs are 050315, 050319, 050401, 050505, 050525A and 050603. 
They are consistent with those of GFM05, apart from two cases. 
For 050401 our L50 measure, where L 50 = L/(10 50 ergs" 1 ), 
is 1405 ± 165, while GFM05 reported 740 ± 100. For 050603, 
we obtained L 50 = 27 06 ± 1470 to be compared with GFM05's 
1200 ± 300. The reason in either case resides in a slightly different 
choice of the 1-s time interval around the peak. GFM05 determined 
this from the 40-350 keV light curve to match the 40-700 keV of 
the BeppoSAX/GRBM, while we used the 15-350 keV. The choices 
of the 1-s time interval turned out to differ by 1-2 s in either case. 
This, combined with the fact that both of these GRBs exhibit a 
sharp peak, turned into the discrepancies provided above. We note 
that in both cases they still lie in the Vr-L region with high Vr and 
high L, consistently with the V/L correlation. 



4 VARIABILITY 
4.1 R01 definition 

The main difference between our data set and those used by R01 
and GFM05 is that our light curves are expressed in background- 
subtracted rates and not in counts. This fact is due to the way BAT, 
which is a coded mask, has been conceived. Hereafter we assumed 
a Poissonian variance for the statistical fluctuations of the light 
curves, as we proved in Appendix A. The formula we used to com- 
pute the variability, according to the R01 definition, is basically the 
same as those of R01 and GFM05, with no extra-Poissonian noise 
term, given that our rates are already background-subtracted. 



V R = 



EN r/\-^N %2 v-^iV 2 2 1 

fel[(L j = i a 'j r i) -£j = i a ij t V j ] 



(2) 



where and bij are the same coefficients as those introduced by 
R01 in their eqs. 6-7. The differences between our formula, eq. 7 
of R01 and eq. 4 of GFM05, are the replacement of the counts Cj 
with the rates rj in the first terms of both numerator and denomina- 
tor, where the original Cj represented the GRB signal, and the re- 
placement of the counts Cj with the statistical noise variances er^ 
affecting the rates Tj in the terms to be subtracted, where the origi- 
nal counts Cj represented the noise. The sum, j = 1, . . . , N, runs 
over the N bins encompassing the GRB time profile. The back- 
ground term Bj in the original formulae of R01 has been set to 
zero, as it has already been removed during the extraction of the 
light curves. 

For each GRB we estimated the smoothing time scale Tj 
(/ = 0.45), defined by R01 as the shortest cumulative time inter- 
val during which a fraction / of the total counts above background 
has been collected. For each GRB we calculated T/ and the cor- 
responding variability Vr, as a function of the binning time. We 
chose the values obtained with the binning time A t that fulfilled 
the requirements reported by GFM05 concerning the ratio At/Tf. 
On one side, when this ratio is too small, the light curve is domi- 
nated by statistical fluctuations, while, on the other side, when the 
binning is too coarse the variability is underestimated. A detailed 
description of these criteria is provided by GFM05. 



4.2 LP06 definition 

Concerning the definition of variability given by LP06, hereafter 
denoted by Vlp, we point out a number of different choices with 
respect to the their analysis. First we estimated Vlp from the 
background- subtracted mask-tagged light curves, while LP06 used 
the raw counts light curves of the 7 Swift /BAT GRBs of their sam- 
ple (Li, private comm.). We assumed no extra-Poissonian variance 
to be subtracted, unlike LP06. We adapted eqs. 1-3 of LP06 ac- 
cordingly and obtained the following: 



Vlp 



Ef =1 [Win 



Vi 



r 2 

' max 



(3) 



(N — 1) r 

where yi is the value for the i-th bin of the reference light curve 
obtained with the Savitzky-Golay filter with a smoothing window 
of T/ (/ = 0.45). W is the same weight as that used by LP06 and 
accounts for the fact that the set of j/i is not completely statistically 
independent from r». As for the determination of the peak count 
rate, r max , we searched the light curve of the same GRB a number 
of times, each time increasing the binning time, until we found the 
peak 5-ct higher than the contiguous bins. This turned out to be 
very accurate, particularly for weak GRBs. In order to comply with 
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the procedure of LP06, N corresponds to the total number of bins 
encompassing the time interval which defines the T90, i.e. from 5% 
to 95% of the total fluence. The values of T90 have been calculated 
using the ftool battblocks. Values of Vlp have been derived from 
the 64-ms light curves. 



5 RESULTS 

Table [T] reports the results of Vr, Vlp, L and T/=o.45 obtained for 
the sample of 41 GRBs. 

5.1 R01 definition 

Significant values of Vr have been obtained for 36 GRBs shown 
in Fig. [T] (circles). In the remaining 5 cases this was not possible 
for different reasons. For GRB 050814 and GRB 050824 we could 
not find any binning matching the requirements mentioned above. 
While for GRB 050126, GRB 050908 and GRB 060512 Vr turned 
out to be consistent with zero within uncertainties. Figure [T] also 
shows the sample of 26 GRBs of GFM05 (squares): the underlu- 
minous GRB 980425, which belongs to the GFM05 sample, is not 
shown because of scale compression reasons; moreover, its uncer- 
tainty on Vr is relatively large. 

We do not show the values GFM05 estimated for six 
SwiftfBAT bursts in common with our sample. Except for the case 
of GRB 050319, our values of Vr for the other 5 GRBs are broadly 
consistent with those of GFM05, some differences being due to a 
different energy band choice (see above). In general, we note that 
our Tf are systematically somewhat higher than those of GFM05: 
this is so because we included low-energy bands, in which GRBs 
are known to last longer. In addition, we know that in some cases 
Vr has a strong dependence on the energy band (GFM05), although 
the definition of Vr by R01 was originally tho ught to account 
for the narrowing of pulses at higher energies iFenimore. et al.l 
ll995MNorris etaljl996h . In the case of GRB 050319 we measured 
Vr = 0.285 ± 0.044, while GFM05 obtained Vr = 0.06 ± 0.03. 
The inconsistency is due to the fact that the original event file, avail- 
able at the time and used by GFM05 to extract the light curve, was 
missing the first sequence of impulses well before the trigger time. 
Therefore, we consider the value reported in this paper as the cor- 
rect one. 

We tested the existence of the V/ L correlation over a number 
of different GRB data sets. Our sample of 36 BAT GRBs shows 
no significant correlation according to Pearson's, Spearman's and 
Kendall's coefficients, whose corresponding no-correlation proba- 
bilities are 72%, 51% and 37%, respectively. However, from Fig.Q] 
we note that in the region of high Vr, and low L, rather unexplored 
by previous data sets (R01; GFM05), there are six GRBs: 050223, 
050416A, 050803, 051016B, 060614, 060729. If one selects the 
BAT GRBs from our sample with L50 > 5, the resulting sample of 
30 GRBs shows a significant improvement of the V/L correlation: 
the probability of no correlation becomes 16%, 5.1% and 3.1%, re- 
spectively. Likewise, if we merge the two samples (GFM05's and 
ours) we obtain similar results: when the 7 bursts with L50 < 5 
are taken out from the total sample of 62 GRBs, the correlation be- 
comes significant with a no-correlation probability of ~ 2 x 10 -4 
according to the non-parametric tests. 

Finally, we calculated Vr in the 25-350 keV energy band, 
i.e. ignoring the lowest energy channel 15-25 keV, of the six low- 
luminosity outliers. The aim was to establish the importance of 
the low-energy channel contribution to the resulting Vr, especially 



when compared with the results of GFM05, whose low-energy 
threshold was 40 keV. We found that in all cases Vr resulted sys- 
tematically higher, although still compatible within uncertainties. 
The only case in which Vr in the 25-350 keV was significantly 
higher than for the whole band was 060614 due to its small statis- 
tical uncertainty. This corroborates the nature of outliers of the six 
GRBs considered: we can rule out that their high values of Vr are 
due to the presence of the low-energy photons not considered by 
previous data sets. 

5.2 LP06 definition 

Significant values of Vlp have been obtained only for 10 GRBs 
shown in Fig. [2] In the remaining 3 1 cases the variability resulted 
consistent with zero within uncertainties (see TableQJ. Despite the 
small number of GRBs with significant Vlp, the correlation ap- 
pears to be significant within 1-2% according to the non-parametric 
tests: 1.1% (Spearman) and 1.6% (Kendall). See Table[2]for further 
details. Figure [2] shows these 10 BAT GRBs as well as the sample 
of 22 GRBs of LP06. Shaded areas show the 1-cr and 2-a regions 
around the best-fitting power law obtained by LP06 using the fitexy 
routine, with a slope of m = 3.25 ± 0.26 and a x 2 /dof = 1.93 
(20 dof). If we ignore GRB 060614, which clearly lies far away 
from any power-law correlation between Vlp and L, and use the 
same routine as LP06, we obtain a best-fitting value of the slope of 
m = 2.3 ± 0.17 and xV dof = 8.5 (7 dof). The X 2 is clearly too 
high and therefore, although the correlation appears to be real, the 
description in terms of a power-law with no sample scatter, as the 
usage of the routine fitexy assumes, is not acceptable. We note that 
this conclusion also holds for the very same result of LP06, whose 
X 2 has a null hypothesis probability of 0.75%. 



6 DISCUSSION 

Interestingly, if one ignores the 6 GRBs from our sample of 
SwiftfBAT with low L, specifically L50 < 5, the remaining ho- 
mogeneous sample of 30 BAT GRBs, for which we could derive a 
reliable estimate of Vr in the 15-350 keV energy band, is fully 
consistent in the Vn-L plot with those from previous detectors, 
thus confirming the existence of the Vr/L correlation. This is re- 
markable, given that BAT is a different kind of 7-ray detector and 
has a different energy band from that of the BeppoSAX/GKBM, 
40-700 keV, whose data mainly comprise the sample of 32 GRBs 
of GFM05. Another important confirmation provided by this BAT 
sample is that the scatter of the correlation originally found by R01 
and GFM05, despite their alternative descriptions of it, is not due 
to the combination of data from different instruments with differ- 
ent effective areas, response functions, statistical noises, and energy 
bands, but it is intrinsic to the correlation. In fact, for the first time 
our data set represents a homogeneous sample of 36 GRBs with 
measured redshift acquired with the very same detector and with 
the very same kind of data for each GRB. 

What is new with this BAT sample is the presence of 6 (out of 
36) low-luminosity GRBs (L 50 < 5). If one ignores GRB 980425, 
a peculiar underluminous and very nearby burst, from the sample 
of GFM05 and R01 it turns out that none of the previous GRBs has 
L50 < 5. This is not surprising, given the unprecedented sensitiv- 
ity of BAT. Therefore these 6 BAT GRBs allow us to test, for the 
first time, whether the correlation holds for low-luminosity GRBs. 
FigureQ]clearly shows that none of them lies where one might have 
expected from the correlation. Instead, they exhibit relatively high 
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Table 1. Variability, according to both definitions considered in the text (Secs. r5~Tl and l5.2t , and peak luminosity for a homogeneous sample of 41 Svvr/f/BAT 
GRBs. 
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a Isotropic-equivalent peak luminosity in 10 50 erg s 1 in the rest-frame 100- 
s" 1 Mpc~\ C m = 0.3, and C A = 0.7. 

values of Vr. This is proven by the correlation coefficients, in par- 
ticular the non-parametric Spearman's r s and Kendall's r, accord- 
ing to which the correlation is significant (5.1% and 3.1% respec- 
tively) or not, depending whether these 6 low-luminosity GRBs are 
excluded or not. This is confirmed by merging our sample of BAT 
with that of GFM05: the correlation is significant, provided that 
low-l uminosity burs t s are e xcluded (see Table[~~). 

iGuidorzi et ail d2007l) have investigated the nature of the 6 
BAT GRBs outliers of the V/L correlation and found strong ev- 
idence that they are also outliers of the anti-correlation, discovered 



1000 keV band, for peak fluxes measured on a 1-s time-scale, Hq = 65 km 



bv lNorris et al.l {2000), between the rest-frame temporal lag and the 
peak luminosity. In particular, they found that these GRBs are char- 
acterised by a small or negligible time lags and a relatively l ow lu- 
minosity. We refer the reader to the paper bv lGuidorzi et al.l d2007t) 
for more details. 

Concerning the definition of variability, Vlp, given by LP06, 
we found that this still correlates with L, although our results dif- 
fer from those by LP06 (see Fig.|2]and Table|2j. In particular, we 
find the description of the correlation in terms of a power law with 
no extrinsic scatter inadequate, given the high values of x 2 /dof 
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Figure 1. Variability Vr vs. peak luminosity L for a sample of 36 long bursts detected by Sw ift/BAT (circles) according to the definition of variability by 
Reicha rt et al] 1200 ll) . For co mparison we show 25 GRBs (squares) from lGuidorzi et all (2005). The shaded areas show the 1- and 1-a regions around the 
best-fit power law obtained by Guidorzi et al. 12006) with the D'Agostini method, with a slope of 1.7. 



Table 2. Correlation Coefficients for different sets of GRBs. 



Set of GRB< a > 






Coefficient (Probability) 








Pearson's r 


Spearman's r s 


Kendall's r 


36 GRBs (Vr vs. 


Lso) 


-0.062 (0.719) 


0.115 (0.506) 


0.105 (0.369) 


30 GRBs (Vr vs. 


L 50 > 5) 


0.261 (0.163) 


0.359 (0.051) 


0.278 (0.031) 


62 GRBsO) (Vr 


vs. L 50 ) 


0.190 (0.139) 


0.315 (0.013) 


0.231 (0.008) 


55 GRBs< b > (Vr 


vs. L50 > 5) 


0.418(1.5 x 10" 


3 ) 0.476 (2.4 x 10" 4 ) 


0.342 (2.3 x 10" 4 ) 


10 GRBs (V LP vs. L 50 ) 


0.536 (0.111) 


0.758 (0.011) 


0.600 (0.016) 



a L 5 = i/(10 50 ergs- 1 ). 

b This sample resulted from the merging of our sample with that of Guidorzi et al. 12005). 



yielded by both samples, ours and LP06's. Regarding our sample 
of 41 BAT GRBs, we find that, unlike the definition of Vr by R01, 
the smoothing filter adopted by LP06 in their definition of VLp cuts 
off the low-frequency variability of GRBs. This results in a selec- 
tion of a smaller sample of GRBs with significant (high-frequency) 
variability: 10 vs. the 36 obtained for the R01 definition. We note 
that GRB 060614 confirms its nature of outlier of the correlation, 
no matter which choice of the definition of variability we adopt 
(Fig.©. 

In general, from Table [2] we note that the Pearson linear cor- 
relation coefficient r is systematically less significant than the 
other two. This supports the finding that the correlation shows a 
clear scatter around the best-fitting power law. Therefore this scat- 
ter must be taken into account properly (e.g. with the D'Agostini 



method), when fitting the data (see D'Agostini 2005 and Guidorzi 
et al. 2006). 



6.1 Low-luminosity GRBs and the Amati correlation 

We tested if the 6 low-lumi nosity GRBs are also outliers of the 
E p ,i-Ei BO dAmati et alj|2002l) (Ei BO is the isotropic energy released 
in the 1 — 10 4 keV rest-frame band) as w ell as of the E p ^-L 
dYonetoku etai]|2004l ; iGhirlanda et alj|2005l) correlations. E 1 = 
E p (l + z) is the intrinsic peak energy of the total spectrum of a 
burst, where E v is the peak of the vF(y) spectrum in the observer 
frame. A correlatio n between temporal variability and E v ,\ was 
originally found by iLlovd-Ronning & Ramirez-R uiz (2002J) for a 
number of bursts with pseudo-redshift derived assuming the vari- 
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Figure 2. Variability Vlp vs. peak luminosity L for a sample of 10 long bursts detected by Swift/B AT (circles) acco rding to the definition of variability by 
iLi & Paczvfiskil fcOOel) . For comparison we show 22 GRBs with significa nt V L p (squares) from|u & PaczyhskJ b006t) . The shaded areas show the 1- and 2-cr 
regions around the best-fit power law obtained by Li & Paczynski 12006), with a slope of 3.25. Solid line shows the best-fitting power law obtained with all of 
the 10 BAT bursts shown here, but GRB 060614; its slope is 2.3 (see Sec.l5~2l 



ability /peak luminosity correlation. Taking into account that E V j 
also correlates with Ei so and with L (isotropic peak luminosity), 
we test whether the breaking of the V/L correlation in the case of 
these 6 bursts is explained by anomalous values of E Pil . 

For two bursts, XRF 0504 16A dSakamoto etal] l2006h and 
GRB 060614 d Amati et al1l2007l) E p ,i has already been reported 
elsewhere. Both GRBs are consistent with the Amati relation. In 
particul ar, XRF 050416A rem arkably confirms it down to the XRFs 
region dSakamoto etai]|2006l) . For the remaining four GRBs, the 
BAT photon spectrum can be fit with a single power law N(E) oc 
E~ Tbat , where That is the photon index. In order to constrain 
E p , we fitted the total spectrum of each burst with a cutoff power 
law by fixing the power law index a to the typical value of 1.0 and 
letting the break energy Eo — E p /(2 — a) free to vary. We took 
the lower/upper limit for Eo from the 90% confidence level inter- 
val on one parameter: if the interval included or lay close to the 
lower (higher) boundary of the BAT passband, we assumed the up- 
per (lower) limit on Eo. Our resul ts are broadly in ag reement with 
the empirical correlation found bv lzhang et alj ( 120071) between E p 
and That- 

Results are reported in Table [3] All of the 6 bursts (or their 
limits) turned out to lie in the 2 sigma region of the Amati relation 
(see Amati 2006). 

We also found that the two bursts with firm E p s as well as 
two with upper limits are consistent with the E Ptl -L correlation, 
while the remaining two lower limits on E Pl i for GRB 050803 
and GRB 060729 are not, as shown in Fig. [3] The better consis- 
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Figure 3. Peak luminosity L vs. rest-frame peak energy j? D ,j of the tota l 
energy spectrum for 5 bursts with firm E Pi i measurements iAmatill2006l) 
and the 6 low-luminosity (L50 < 5) GRBs (empty circles) of our Swift/BAT 
sample. 



tency with the E Pi i-Ei BO than with the E Ptl -L correlation can be 
explained with the fact that E Pl \ better correlates with the time in- 
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Table 3. Intrinsic peak energy E p ; of the total spectrum for the subset of 6 
low-luminosity GRBs of our sample. Tbat is the photon index of the total 
photon spectrum (N(E) oc E~ f bat) when this is fit with a single power 
law in the BAT energy band. Limits are given at 90% confidence level. 



GRB 



Tbat 



B p ,i (keV) E iso (10 52 erg) 



050223 

050416A< b ' c ) 

050803 

051016B 

060614( d ) 

060729 



1.90±0.16 (a) 

1.58 ± 0.09 
2.13 ± 0.27 

1.62 ±0.18 



< 127 
25.1 ±4.2 

> 103 

< 125 
55 ±45 

> 80 



0.12 ±0.02 
0.12 ±0.02 
0.20 ±0.03 
0.14 ±0.04 
0.25 ±0.10 
0.27 ±0.05 



a In agreement with lPage et alj i2005h . 

b from lAmatil J2006l) . 

c from Sa kamoto et all l200r|) . 

d from lAmati et alj 120071) . 
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Figure 4. Rest-frame duration Tgo.rcst vs. peak luminosity for all the 
41 Swift/BAT bursts reported in Table [T] Empty circles show the 6 low- 
luminosity (L50 < 5) high-variability GRBs. 



7 CONCLUSIONS 

We tested the variability/peak luminosity (V/L) correlation with 
a homogeneous sample of 36 GRBs detected with Swift/BAT in 
the 15-350 keV energy band with firm redshift. We adopted two 
different definitions of variability: that by Reichart et al. (2001; Vr) 
and that by Li & Paczynski (2006; Vlp), which differs from the 
former for a different smoothing filter. The most interesting results 
have been derived with Vr. If we consider only the GRBs with peak 
luminosity L comparable with those of previous samples, i.e. L > 
5 x 10 50 erg s _1 in the rest-frame 100-1000 keV energy band, we 
confirm the correlation and i ts intrinsic dispersion around the best- 
fitting power law obtained bylGuMOTzretayl^OO^): m = 1.7±0.4 
(L oc V m ) and <r logL = 0.58±°;J|. 

Interestingly, all of the 6 low-luminosity GRBs detected by 
Swift/BAT turn out to be outliers to the V/L correlation, show- 
ing higher values of Vr, than expected. This does not contradict 
the results from previous samples of GRBs with known redshift. 
Instead, we are led to conclude that the correlation does not hold 
any more for low-luminosity GRBs. We also fi nd that these bursts 
are consistent with the E p a-E\ so correlation jAmati etaL 20021) 



and four of them also with t he E Ptl -L correlation 1 Yonetoku et al.1 

l2004l ; lGhirlanda et alj2005l) . 

Unlike the results obtained bv lLi & P aczynski (2006), we do 
not find evidence for a tighter correlation using Vlp instead of 
Vr. Rather, fewer GRBs appear to have a significant measure of 
Vl p; we ascribe th i s to th e fact that the smoothing filter adopted 
by iLi & PaczvriskH (2006) to construct the reference light curve 
with respect to which the variability is computed, only selects high- 
frequency variability. 
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tegrated released energy, as proven also by the sc atter of the corre- 
lation between L and i5i so (Ghirlanda et al. 2005). 

We also tested whether the duration of these events correlates 
with their peak luminosity. To this aim, in Fig. [4] the rest-frame 
T9o,rest = Tgo/(l + z) is plotted against L50 for the entire sample 
of 41 Swift/BAT GRBs considered. X90 is the time interval col- 
lecting from 5% to 95% of the total fluence in the observer frame. 
For each burst we used the value published by the BAT team in 
the refined GCN circulars. Empty circles correspond to the 6 low- 
luminosity GRBs with a significant measure of variability. Appar- 
ently there is no hint for correlation and also no evidence for a 
different behaviour of the 6 low-luminosity GRBs with respect to 
the others. The result does not change in essence when we replace 
Tgo.reat withTgo. 

We conclude that the fact that the variability of these 6 low- 
luminosity high- variability GRBs does not correlate with the peak 
luminosity is not connected with their E Pl \, which correlat es with 
E'iao as almost all of the long GRBs with known redshift lAmatil 
120061) . 



APPENDIX A: STATISTICAL NOISE OF BAT 
MASK-TAGGED LIGHT CURVES 

We report the analysis performed on the BAT mask-tagged light 
curves of the GRBs considered in this work, aimed at studying the 
statistical noise. As the GRB itself is characterised by intrinsic tem- 
poral variance which is unknown a priori, we limited to the pre- and 
post-burst regions of the light curves, where the background is sup- 
posed to be the dominant source of statistical fluctuations. In order 
to make sure that we excluded the entire light curve of the GRB, 
we binned it spanning very different integration times (from 64 ms 
to 32 s) and checked that no trend in the residuals was visible. 

Let n and a Tj be the count rate and its uncertainty, respec- 
tively, of the i-th bin of a 64-ms mask-tagged BAT light curve. This 
light curve is the result of the BAT standard pipeline already sum- 
marised in Sec. [2] (see also Barthelmy et al. 2005). Uncertainties 
oy. (i = 1, . . . , TV, where N is the total number of bins of the 
selected portion of light curve) are calculated by propagation of er- 
rors, starting from the raw counts assumed to be affected by purely 
Poissonian noise through the ftool batbinevt. 
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We tested the following null hypothesis: each n is a single re- 
alisation of a normal random variable with null expected value and 
sigma equal to a Ti : N(0, cr ri ). Little can be inferred on a random 
variable from a single realisation. However, as long as this hypoth- 
esis is true, the various ri/a ri (i = 1, . . . , N) can be seen as dif- 
ferent realisations of the same random variable, r n , hereafter called 
"normalized rate", which has a standard normal density: N(0, 1). 

So we studied the observed distribution of r n for each single 
light curve removed of the GRB profile. We fitted the observed 
distribution with a Gaussian N(fi, a). 

In particular, we are interested in constraining the possible 
presence of any additional source of statistical noise (e.g. instru- 
mental) to the Poissonian one. 

More generally, should the various fi fluctuate more than a ri , 
so that the true variance is (1 + / np ) a 2 . , where / np is the fraction 
of additional non-Poissonian variance, the resulting a should be 
greater than unity. More precisely, we should find a 2 = (1 + / np ). 

Therefore we fitted the observed distribution of r n with 
N(fj,, a), first by imposing a — 1. In every case we found accept- 
able x 2 values, confirming that no evidence for additional noise has 
been found. 

In particular, we were interested in setting a limit to / np with a 
given confidence level. Following Papoulis & Pillai (2002; p. 313— 
314), in the case of unknown fx we used the sample variance s 2 
defined as: 



N — 1 



(Al) 



where r n< i = Tija ri is the single realisation of r n and r n is the 
mean value. The random var iable (TV - l)s 2 /a 2 follows a % 2 {N - 
1) distribution, so that we can constrain a 2 , i.e. (1 + / np ), through 
the following: 



l + /n 



a 2 < 



Ds 2 



X 2 S/2 (N- 



(A2) 



at (1 — S) confidence level; xt( n ) is me u percentile of the X 2 ( n ) 
distribution. In most cases N was big enough (> 10 3 )to ensure the 
following approximation: 



/np < 



l + 2l- 



5/2 



(N-l) 



- 1 



(A3) 



where z u is the u percentile of the standard normal density. 

We show the example of GRB 050401. The distribution of 
r u ,i, N = 7065, can be fit with a Gaussian N(0, 1), \ 2 /Aoi = 
32.2/28, as shown in Fig.lATl 

The sample variance resulted s 2 = 0.999. The consequent 
upper limit on f np turns out to be 2.7% (4.3%) at 90% (99%) con- 
fidence level. 
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